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ABSTRACT 

Secondary anisotropies of the cosmic microwave background (CMB) can be detected by using 
the cross-correlation between the large-scale structure (LSS) and the CMB temperature fluc- 
tuations. In such studies, chance correlations of primordial CMB fluctuations with the LSS 
are the main source of uncertainty. We present a method for reducing this noise by exploit- 
ing information contained in the polarisation of CMB photons. The method is described in 
general terms and then applied to our recently proposed optimal method for measuring the 
integrated Sachs-Wolfe (ISW) effect. We obtain an expected signal-to-noise ratio of up to 8.5. 
This corresponds to an enhancement of the signal-to-noise by 23 per cent as compared to the 
standard method for ISW detection, and by 16 per cent w.r.t. our recently proposed method, 
both for the best-case scenario of having perfect (noiseless) CMB and LSS data. 
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1 INTRODUCTION 

The low-redshift large-scale structure (LSS) changes the cos- 
mic microwave background (CMB) fluctuations in various ways. 
Such secondary effects on the CMB are, for exampl e, the inte- 
grated Sachs- Wolfe ( ISW) effect dSachs & Wolfdl 19671) . the Rees- 
Sciama (RS) effect JRees & Sciamalll968l) , gravitational lensing 
( [Lewis &~C hallinor 2 006h. and the Sunyaev-Zel'dovich (SZ) effect 
■ ISunvaev & Zeldovichlll972lll98Cl) . By studying these signals, we 
can obtain valuable information about our Universe. The ISW ef- 
fect, for example, provides independent evidence for the existence 
of dark energy. Unfortunately, unless the spectral signatures of the 
signal differ from the ones of the primordial CMB, it is difficult 
to detect them. The reason is that the primordial CMB fluctuations 
created at the time of last scattering are much stronger than the 
secondary temperature anisotropies. The usual method for detect- 
ing secondary anisotropies in the CMB is via cross-correlating the 
CMB temperature maps with LSS data such as the galaxy density 
contrast. Since secondary anisotropies in the CMB are created by 
the LSS, there is a significant cross-correlation between the two. 
In contrast, the primordial CMB fluctuations should not be corre- 
lated with the LSS. By performing the cross-correlation analysis, 
one can therefore separate signatures of the presence of these sec- 
ondary anisotropies from the primordial fluctuations. 

In the standard cross-correlation method, first described by 
Boug hn et al. Ul998l) . the observed cross-correlation between LSS 
and CMB data is compared to its theoretical prediction. This 
method has been extensively used to detect the ISW effect. Some of 
the most recent studies are bv lHo et all ( l2008h . lGiarinantonio et al.l 



d2008h . iRassat et alj ( l2006h . and IBoughn & Crittendenl j2004h . 
Since the theoretical cross-correlation function is by construction 
an ensemble average over all possible universes, fluctuations asso- 
ciated with the specific realisation of the LSS in the observed Uni- 
verse act as a source of noise in the detected signal in the standard 
method. 

In IFrommert et al. I fe008l) . we suggested a method for re- 
ducing this source of uncertainty, which we will refer to as the 
optimal temperature-only method. Instead of comparing the ob- 
served cross-correlation function with its theoretical prediction, we 
use an optimal matched filter in order to detect an ISW template 
in CMB data . Similar schemes were independ ently proposed by 
IZhand J2006I) . iHernandez-Monteagudd d2008ri and iGranett et afl 
(2008). Optimal matched filters have also been used to study 
other secondary effects on the CMB. The first of these studies 
expl ored the detectability of the kinetic SZ effect of galaxy clus- 
ters (Haehnelt & Tegmark T996h . later works on the kinetic SZ ef- 
fect and the RS effe ct are for example by ISchafer et alj < 2006). 
IMaturi et al] fe007al) . iMaturi et al] J2007bh . and IWaelkens et al] 
j2008l) . 

However, in both the standard method and the optimal 
temperature-only method, the main source of uncertainty in the de- 
tection of the secondary signal comes from chance correlations of 
primordial CMB fluctuations with the LSS. In this work, we present 
a method which exploits polarisation information in order to reduce 
not only the noise from the specific LSS realisation, but also the 
noise coming from primordial CMB temperature fluctuations. This 
method can be applied generically to the detection of all secondary 
effects. It is based on the fact that the polarisation measured in the 
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CMB contains information about the primordial temperature fluc- 
tuations. We use the observed E-mode polarisation map, which we 
translate into a temperature map using the TE cross-power spec- 
trum. The obtained temperature map is then subtracted from the 
observed temperature map, and hence no longer contributes to the 
noise budget of the detected signal. Once an E-map has been mea- 
sured to a good accuracy, this will significantly enhance the signal- 
to-noise ratio of the detection of secondary effects. The first all- 
sky measurement of polarisation with high fi delity is expe cted to 
be provided by the Planck Surveyor satellite (Tauber 200G), to be 
launched in 2009. 

Our optimal polarisation method builds on the optimal scheme 
to detect LSS signatu res in CMB data, which we developed in 
iFrommert et al.l d2008h specifically for the ISW effect. Note that 
this method assumes a Gaussian data model, hence it is very well 
suited for the ISW effect, whereas one might need to extend it into 
the non-Gaussian regime for other effects, such as the RS effect, 
the kinetic SZ e ffect or lensing. Th is can be done using informa- 
tion field theory (Enfil in et alj2008h . but is beyond the scope of this 
work. Here we show how to use the information contained in po- 
larisation data within the framework of a Gaussian data model and 
leave the extension to more complicated models for future work. 

When applying our method to ISW detection, we obtain an 
expected signal-to-noise ratio of up to 8.5. This corresponds to an 
enhancement of the signal-to-noise ratio by 16 per cent w.r.t. the 
optimal temperature-only method, independent of the depth of the 
galaxy survey considered. In comparison to the standard method, 
the signal-to-noise ratio is enhanced by 23 per cent for a full-sky 
LSS survey that goes out to redshift 2. Both of these comparisons 
have been made for the best-case scenario of having perfect (noise- 
less) CMB and LSS data. 

Using polarisation data to reduce the noise in the detection of 
secondary effects was first proposed b y Robert Critten den, follow- 
ing a suggestion from Lyman Page {Crittenden 2006). He already 
derived the reduced temperature power spectrum, which we show 
in Figure Q] and roughly estimates the improvement of the signal- 
to-noise ratio for ISW detction to be around 20 per cent, which we 
confirm with our calculations. 

Our article is organised as follows. In sect ion [2] we describe 
the optimal method derived in IFrommert et alj (|2008|) in general 
terms. In section [3] we then show how we can reduce the noise 
coming from primordial temperature fluctuations by using polar- 
isation data. In section [4] we apply the method to the ISW effect. 
We conclude in section[5] 



2 OPTIMAL METHOD FOR THE DETECTION OF 
SECONDARY EFFECTS ON THE CMB 

In IFrommert et alj J2008), we derived an optimal method for the 
detection of secondary temperature anisotropies in the CMB using 
as example the ISW effect. In this section we briefly review this 
method, which we refer to as the optimal temperature-only method. 

Let's assume that we know the LSS well enough to create a 
template TV of the secondary signal T s that we would like to detect 
in the temperature fluctuations, for example the ISW signal Ti sw . 
Here, Tx with any index X denotes the function Tx '■ S 2 — > R, 
which we regard also as an element of a function- vector space. The 
data d we measure are the observed CMB temperature fluctuations 
T b s . Our data model is then 

d = T bs 



Symbol 


Definition 


^cmb' ^cmb 


cosmological CMB temperature and polarisation 


T s , E s 


real secondary signal that we are trying to detect 


T T ,E T 


signal templates for temperature and polarisation 


T fg , E fg 


residual galactic foregrounds alter foreground removal 


Tdet > Edet 


detector noise 


Tobs > -^obs 


(Tcmb + Tfg + Tdet), {Ecmh + -Efg + -^dot) 


AT oba , A_E obs 


(T bs — TV), (E obs — E T ) 


T- 

+ ISW 


fluctuations created by ISW effect 


T ■ 

± prim 


(T cm b Tj sw ) 



Table 1. Summary of defined symbols 

+ Tf g + Tdet 
= T T + (T cm b — T T ) + Tfg + Tdet 

= T T + AT obs , (1) 

where T cm b denotes the cosmological CMB temperature fluc- 
tuations, Tfg are residual galactic foregrounds after foreground- 
removal, and Tdot denotes the detector noise. Note that (T cm b — 
TV) = (T cmb - T s ) + (T s - T T ) contains the the CMB fluctuations 
other than the signal we are after, (T cm b — T a ), and the uncertainty 
in the template w.r.t. the signal, (T s — TV), coming from our igno- 
rance of the full distribution of the matter in the Universe. Note that 
for simplifying the notation, we have redefined T = (T — To) /To, 
where To denotes the monopole of the CMB temperature fluctua- 
tions. An overview over the above definitions can be found in ta- 
ble [j] 

We now approximate the distribution of AT b s by a Gaussian 
around zero. That is, we write the probability density function of 
T bs given the signal template TV and the cosmological parameters 
p, the likelihood, as 

T(T ob s | T T , p) = Q (T bs - T T , C) . (2) 
Here we have defined 

^C),-^= r {-^C- X ) (3) 

to denote the probability density function of a Gaussian distributed 
vector x with zero mean, given the cosmological parameters p and 
the covariance matrix C = (xx^)> where the averages are taken 
over the Gaussian distribution Q{\,C). Note that in general the 
covariance matrix depends on the cosmological parameters, which 
is not explicitly stated in our notation. A daggered vector or matrix 
denotes its transposed and complex conjugated version, as usual. 
Hence, given two vectors a and b, a must be read as the tensor 
product, whereas a* b denotes the scalar product. Note that in eq. 
([2} the signal template TV may depend on the cosmological param- 
eters p as well. 

Let us briefly address the question of how to create the tem- 
plate T T . When writing down the likelihood in eq. l[2}, we have 
implicitely assumed that the template TV is the mean of T b s w.r.t. 
the probability distribution given in eq. {2j. This probability distri- 
bution is conditional on the template TV, or, in other words, condi- 
tional on the LSS data 5 g , from which we have created our template 
according to some prescription. Note that usually <5 g denotes the 
galaxy density contrast, but we use it to denote the LSS data in a 
more general sense here, which could also be lensing information, 
for example. 

In the following, we assume that the signal T s = R8 m is 
given by a linear operator R applied to the matter density contrast 
Sm- For the ISW effect, the operator R is explicitely derived in 
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(4) 



From mert et al. We can then write 

T T = {T' b s )p(T obs | s g ,p) 

~ (r s )p(T s | Sg,p) + ({Tcmb ~ r s )) P (( Tcmb _ Ts) | 5gjp ) 

+ ( T fg>P(T f g | i s ,p) + {2dct)p(T dct ) 
= R(Sm)p(S m | S S ,J>)> 

where we have used that T s , (Tcmb — T s ), Tf g and Td e t are approxi- 
mately stochastically independent in the first step, and that the three 
errors have vanishing means, ((T cmb - r s ))p((T cmb -T a ) | s g , P ) = 
(Tig}p(T ts \6 g , P ) = <Tdet)p(T dct ) = 0, in the second step. 
For the ISW effect, (T cmb - T s ) = (T cmb - T isw ) = T prim 
are simply the primord ial fluctuations, which do have zero mean 
dFrommerte t al. 2008). For other secondary effects, ((T cm b — 
T s )}p((T cmb -T a ) | s„,p) = is probably still a reasonably good ap- 
proximation. In the last step, we have pulled the operator R out of 
the mean. 

We see that for creating the signal template TV, we need the 
mean of the matter density contrast conditional on the LSS data, 
(Sm) p(s m \s , P )- In the simplest case of having a Gaussian likeli- 
hood and Gaussian prior for 8 m , this is given by the Wiener filter. 
Again this is a very good approximation for the ISW effect, which 
is present on very large scales, on which structure growth is still 
linear. For other effects such as the kinetic SZ effect, the RS effect 
or lensing, the Gaussian approximation for S m may not be very 
good (thus also the Gaussian approximation for AT b s may not be 
good), and one would have to consi der non-Gaussian data mod- 
els using information field theory ( lEnfilinetal1l2008l) . However, 
in this work we will use the Gaussian data model and leave ex- 
tensions to non-Gaussian models for future work. Note that, when 
choosing the template as in eq. {4]l, the the latter is uncorrelated 
with AT b s = (T bs — TV) (w.r.t. the probability distribution in eq. 
(f2]l), as can be easily shown. 

In order to see how well we can recover such a signal template 
from the CMB data, we put an amplitude A T in front of the signal 
T T in eq. {2), and try to estimate its value from the data (the true 
value of this amplitude is one, of course): 

P(T obs I A T ,T T ,p) = G(T ohs - A T T T , C). (5) 
The maximum likelihood estimator A T for the amplitude A T is 

£;(2Z + 1)^- 



E ; (2Z + 1)-^ 



(6) 



In the second equality, we have assumed that the knowledge of the 
secondary anisotropy template is equally good in any direction, so 
that the template uncertainty matrix is isotropic and fully described 
by its spherical harmonics power spectrum. We will use this as- 
sumption also in the following. This permits us to evaluate the ex- 
pressions in spherical harmonics space in the second step. We have 
used the following definitions of the power spectra and their esti- 
mators (we use a hat to denote estimators) 



f~iX,Y I X Y*\ 

— \ a lm a lm ll 



Cf 



cf 



c 



21 



— ^2 Re (ahnafm^J , 



Cf 



Cf 



(7) 
(8) 

(9) 
(10) 



where the a; m are defined by an expansion into spherical harmonics 

Y lm : 



af m = / mT x {n)YC m (n). 



(11) 



The power spectrum cf* T ° bs denotes the spherical harmonics space 
version of the covariance matrix C . We calculate the variance of the 
amplitude estimator to be 



a A = ([A T - (Ar)cond ) ) 



(tSc^Tt) 1 = (]T(2« + 1) 



Cf 



C AT obs 



(12) 



where we have again evaluated the expressions in spherical har- 
monics s pace in the last step, an d we have used the notation intro- 
duced in Frommert et al. where the index "cond" indicates 
that the average is taken conditional on the signal template T T , i.e. 
over the probability distribution given in eq. {2j. We can now define 
the signal-to-noise ratio as follows 



Cf 



(13) 



where the index t indicates that this is the signal-to-noise ratio one 
obtains for the optimal temperature-only method. This signal-to- 
noise ratio depends on the actual realisation of t he matter distri- 
bution in our Universe via the estimator C; Tt . In Fromm ert et al.l 
(2008), we showed that for the ISW effect we obtain on average a 
signal-to-noise ratio of about 7, if we assume an ideal LSS survey 
which covers the whole sky and goes out to a redshift of about 2. In 
comparison to the standard method, this is an enhancement of the 
signal-to-noise ratio by about 7 per cent. 



3 REDUCTION OF THE PRIMORDIAL NOISE USING 
POLARISATION INFORMATION 



With the method suggested in lFrommert et alJd2008l) . we were able 
to reduce the low redshift cosmic variance effect in amplitude esti- 
mates of secondary signals, i.e. we reduced the noise coming from 
the specific realisation of LSS in our Universe. Now we show how 
even the noise coming from primordial temperature fluctuations 
can be reduced. The idea is that since the temperature and polarisa- 
tion maps of the CMB are correlated, the polarisation contains in- 
formation about the temperature fluctuations. After extracting this 
information from the polarisation data we know a part of the tem- 
perature map, which we can remove from the data before trying 
to detect the signal. In other words, we make our amplitude esti- 
mate of the secondary signal conditional on the known part of the 
temperature fluctuations. 

To include the information contained in the polarisation data, 
we enlarge our data vector d to include the observed E-mode polar- 
isation map -Bobs as well: 



d = (Tobs, Sobs) , 

or, in spherical harmonics space 



d 
a lrr, 



r o b 3 ^obs^ 
!m > u lrn I 



(14) 



(15) 



Note that with the map E a ^ s , we are referring again to the abstract 
element of a function-vector space space, which contains all the 
information on the observed E-mode. When evaluating the abstract 
expressions obtained in the following, we use the representation 
of iJ bs in spherical harmonics space, consisting of all coefficients 
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In principle, it is possible that the secondary effect we are 
looking for is also present as a small signal in the polarisation 
data. If the temperature anisotropies created by the secondary ef- 
fect exhibit a quadrupole component at the time of reionization, 
this quadrupole wil l be rescattered by free electrons and create a 
polarisation signal I Zalda rriagal l 1997). How ever, for the ISW this 
effect has been proven to be small dCoorav & Melchiorrill2006l) . It 
should also be small for the RS effect, lensing and the kinetic SZ 
effect, the highest contributions of which are on relatively small 
scales. Thus, as a first approximation we assume that the polarisa- 
tion data do not carry any signal of the effect we want to detect. 
Our signal template r is then 



a ln 



(TV,o) f , 



and the data model becomes 



d = 



Tobs 
E bs 



T T + Ar o b s 

-Eobs 



(16) 



(17) 



The observed E-map, E a ^ s — _E cm b + Ef g + E dct , consists of 
the cosmological E-mode fluctuations i5 cm b, residual galactic fore- 
grounds after foreground removal E{ g , and the detector noise E de t- 
Assuming again Gaussianity, we can write down the likelihood 



p(d\T, P ) = g{d-r, c), 

where the covariance matrix C is 

C={(d-T)(d-T)% oad , 



(18) 



(19) 



and we have redefined the index 'cond' to denote the average 
over the probability distribution in eq. J 1 8 fc . In spherical harmonics 
space, the covariance matrix C is block-diagonal with the blocks 
being 



C(0 = ( c AT obe ,E 



AT oba ,B obs 
C E obB 



Therefore, the likelihood factorises: 
P(d\r,p) = Y]_G(af m - aJ m ,C(l)). 

l,m 



(20) 



(21) 



the likelihood of the reduced temperature map 



P(T led \T T ,p) 



TV ^rcds 



(24) 



Note that it is straightforward to derive the factorised likeli- 
hood also for the case that we do have a non-zero signal template 
E T for the polarisation part. In that case, the covariance matrix C(l) 
is slightly changed, as well as the definitions of the reduced tem- 
perature map and power spectrum, and we can no longer neglect 
the polarisation part of the likelihood. Please refer to Appendix lAl 
for details. 

Let us pause for a second and have a closer look at the quan- 
tities defined in eq. d22t . What we have effectively done is the fol- 
lowing. We have a polarisation map af^ bs , which is correlated with 



the temperature fluctuations a,. 



lm 

via C, 



. That is, the 



polarisation map contains information about the temperature map, 
which we can translate into a 'known' part of the temperature map 
using the prescription ^cf^ T ° bs ' E ° b!S /Cf oba ^j a E ° ba - This known 
part of the temperature map is subtracted from the observed one, 
and we work only with the remaining unknown temperature fluctu- 
ations in which we try to detect our signal template. 

The reduced temperature map fluctuates around our signal 
template T T only with the variance C[ cd , which is smaller than 
the full variance C* ; AT ° ba of the observed temperature map. This 
reduced variance is the uncertainty going into our signal detection 
problem now, rather than the full variance of the original tempera- 
ture fluctuations. 

In order to see this, let us again put an amplitude in front of 
the signal template in eq. d24b . and estimate it from the data using 
a maximum likelihood estimator: 



T Icd 'C Icd T T _ t / , 

T ^c; c )r T £,(2Z + i)^ 



(25) 



Here, the last expression is in spherical harmonics space. The vari- 
ance of A T is now 



When inserting the inverse of the covariance matrix C(l), it is pos- 
sible to rewrite the likelihood as a product of a reduced temperature 
part and a polarisation part. To this end, let us define the reduced 
temperature map and power spectrum 



C7 



Cf 



lm ' 
2 



ci 



With these definitions, the likelihood becomes 



P(d\T,p) = H[g(a 



I m 



1,771 



(22) 



(23) 



as we prove in Appendix lAl Now our goal is to find the signal tem- 
plate T T in the CMB data. The polarisation part of the above like- 
lihood, Q{af^" , Cf° bB ), does not depend on the signal template, 
nor does the reduced temperature part explicitely depend on -Eobs- 
In other words, the observed E-map does not contain relevant in- 
formation any more after introducing the reduced temperature fluc- 
tuations. Thus, we can marginalize over it, and continue only with 



c 



' £Ted 



(26) 



and hence the signal-to-noise ratio becomes 



TV 



pol 



E 



(2l + l)Cf 



(27) 



Note that we have added the index "pol" to indicate that this is the 
signal-to-noise ratio one obtains when using the polarisation data to 
reduce the variance. Comparing the signal-to-noise ratio in eq. Ml\ 
with the one in eq. d 13b . we see that by including the information 
contained in the polarisation data, we reduce the variance in every 

mode by the term (cf T ° bB ' B ° ba ) /Cf oba . 

Let us now get an impression of how much the variance gets 
reduced for the different multipoles. To this end, we neglect the 
detector noise Tdet and Edct, and the foreground noise Tf g and 
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Figure 1. Reduction of the variance in the detection of secondary 
temperature signals by using the information contained in polarisation 
data. Shown are the CMB temperature power spectrum C ; Tcmb (solid), 
and the template-free part of the reduced temperature power spec- 



trum C 



C, 



T cmb ,-Ec„ 



/C ; cmb (dashed), together with the 
part of the CMB power spectrum coming from the 'known' part of 
the temperature fluctuations which we infer from the polarisation map, 



C 



/C ( T cmb (dotted). 



Ef JE which allows us write 



CT 



O 



c 



2C, 



c; 



E,_ii 



T T ,E C „ 
I 



(28) 
(29) 

For the ISW ef- 



We furthermore neglect the cross-term C ( 
feet, we have verified numerically that it is negligible. For the ki- 
netic SZ and RS effects, the template itself is so small that we 
can also certainly neglect (jf T ' EcTab . Then, the reduced tempera- 
ture power spectrum defined in eq. d22t becomes 



or 



o 



2c; 



T CIah ,T T 



c 



(30) 



In Fig.[T] we plot the template-free part of the reduced temperature 



C, 



/c, 



power spectrum C l 

have not included the template-dependent terms 



(note that we 



-2C7 cmb ' JT and 

Cj T in the plot), which gives us an impression of how the variance 
coming from primordial temperature fluctuations is being reduced 
by including polarisation data. The variance will be further reduced 
by working conditional on the signal template TV, which is en- 
coded in the terms — 2C ; Tcmb ' Tt and Cj T , and already described in 



From mert et al. | j2008h . We also plot the original CMB power spec- 
trum (7 ; Tcmb and the difference to the reduced one for comparison. 
We have assumed a flat ACDM model with the parameter values 
given bv lKomatsu et al.l J2008l) , table 1 (Q. b h 2 = 0.02265, fi A = 
0.721, h = 0.701, n s = 0.96, r = 0. 084, q 8 = 0.817), and 
used CMB EASY (www . cmbeasy . org , lDoranl2005l) for obtain- 
ing the respective spectra. 

In Fig. [2] we plot a realisation of the original temperature map 
Tcmb (top panel), the reduced temperature map T rc d (middle panel) 



In reality, galactic E-mode foregrounds E[ g are likely to be the limit- 
ing factor in the improvement of the detection significance coming from 
including polarisation data. We comment on this at the end of this section. 




Figure 2. Realisation of the original CMB temperature map T cul i, (top 
panel), the reduced temperature map T rcd (middle panel) and the difference 
between the two for comparison (bottom panel) in fiK. We have chosen the 
same colour range from — 500(j,K to 500fiK for all maps. 



and the difference of the two, lcf ,T ° be ' E ° ha /Cf° j u, lm 
comparison (botto m panel). The realisa tions were created using the 
HEALPix package dGorski et al]2005h . 

Note that all of what we have done works equally well for 
reducing the E-mode polarisation map when trying to detect a sec- 
ondary signal contained in the polarisation data. One has to simply 
exchange the roles of T an d E in the derivation. This was partly al- 
ready done bv lJaffd d2003h . who used the information contained in 
the CMB temperature map for predicting a polarisation map from 
it. The equivalent plot to Fig. Q] for this scenario is given in Fig. 
[3] The likelihood for the case of simultaneously detecting a tem- 
perature template T T and a polarisation template E T is derived in 
Appendix [A] 

In practice, the accuracy to which we can measure the E-map 
is limited by galactic foregrounds Ef g , the most important of which 
are synchrotron radiation and dust emission of the Milky Way. Un- 
certainty in the measured E-map makes the reduction of the temper- 
ature power spectrum less efficient, because the power contained in 
the foreground noise, C\ ts , enhances the observed E-mode power 



for 



spectrum C t 



Qb +C , f fE +C, Bdet . 



The prediction of a re- 



alistic signal-to-noise ratio for our method would require a detailed 
study of foreground effects, detector noise, and scanning strategies, 
which is beyond the scope of this work. 
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Figure 3. Reduction of the variance in the detection of secondary polar- 
isation signals by using the information contained in temperature data. 
Shown are the CMB E-mode power spectrum C i B ° mb (solid), and the 
template-free part of the reduced E-mode power spectrum C ; Ecmb _ 

(dashed), together with the part of the CMB 



C 



vcf- 



power spectrum coming from the 'known' part of the E-mode fluctuations 

(dot- 



■*) /C 



which we infer from the temperature map, \^C l ' 
ted). 

4 EXAMPLE: THE ISW EFFECT 



Let us now apply our method to the ISW effect. That is, our sig- 
nal template T T is now an ISW template which we obtain from a 
Wiener filter reconstruction of the LSS, wh ich can be shown to b e 
optimal for the purpose of ISW detection jFrommert et alj |2008). 
We assume the best-case scenario of having perfect (noiseless) LSS 
and CMB data. In other words, we neglect the detector noise Tdet 
and -Bdct, which is safe on the largest scales, where cosmic vari- 
ance dominates (Afshordi 2004). We furthermore neglect residual 
galactic foregrounds Tf g and Ef & as well as the shot-noise in the ob- 
served galaxy distribution, and assume that we have an ideal galaxy 
survey that covers the whole sky and goes out to a redshift of at 
least two. Then our signal template is exact, T T = T s = Ti sw , 



and the residual (T cm b — Tis 



is simply given by the 



primordial CMB fluctuations, which are created at the surface of 
last scattering (we have ignored other secondary effects here). We 
further assume T sw to be uncorrelated with the primordial fluctu- 
ations Tp r i m , which is a safe assumption bec ause they are created 
on very different scales (Boug hrTet al.l 1 1998h . We can then write 



C 



c 



= c, 



The signal-to-noise ratio for the detection of the ISW signal, 
eq. < !27b . then reduces to 



Nj 



pol 



(2l + l)C? 



i c; 



c 



(3D 



As we said before, the signal-to-noise ratio depends on the specific 
LSS realisation in our Universe via C ; T ' aw . We can infer its prob- 
ability distribution from the distribution of Ti sw by using the cen- 
tral limit theorem for the distribution of ( S/N) and deriving t he 
distribution for S/N from that (see also iFrommert et al.l l2008pl 
We then average the signal-to-noise ratio over this probability dis- 
tribution in order to compare it to the signal-to-noise ratio of the 
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Figure 4. Comparison of the cumulative signal-to-noise ratios for z max = 
2. Top panel: Average signal-to-noise ratio of the optimal polarisa- 
tion method (S/N)^^ (solid), of the optimal temperature-only method 
(S/N)t v (dashed), and signal-to-noise ratio of the standard method 
(S/N) B t (dotted) versus the maximal multipole considered in the analy- 
sis. Bottom panel: Ratio of the signal-to-noise of the optimal polarisation 
method with the one of the standard method (solid) and with the one of the 
optimal temperature-only method (dashed). 



standard method and the average signal-to-noise ratio of the op- 
timal temperature-only method, both described in IFrommert et al.l 
(2008). Recall that the signal-to-noise ratio one obtains for the stan- 
dard method is given by 



= E 



(2* + l)C, J 



c 



(32) 



The cumulative signal-to-noise ratios versus the maximal multipole 
'max used in the analysis are plotted in Fig. [4] Here we have as- 
sumed the ideal galaxy survey described above. We see that includ- 
ing the polarisation data in the analysis increases the signal-to-noise 
ratio by 16 per cent as compared to the optimal temperature-only 
method, and by 23 per cent as compared to the standard method. 
Note that we only included the linear ISW effect in Fig. [4] Beyond 
a multip ole of about I « 100, non-linear effects start to play a cru- 
cial role ( ICooravll2002l) . which could change the plot for I > 100. 
However, we see that for the linear ISW effect, there is hardly any 
contribution for such high multipoles. 

Let us now look at the enhancement of the signal-to-noise ra- 
tio for shallower LSS s urveys. We use the same approximation as in 
IFrommert et al.l J2008), i.e. we introduce a sharp cut-off in redshift 
and redefine everything beyond that redshift as primordial fluctua- 
tions. This introduces a correlation between what we consider the 
ISW and primordial fluctuations, which we would not have if we 
had used a proper Wiener filter based template T T for redefining 
Ti sw . However, for getting a rough picture of the redshift depen- 
dence, this approximation is good enough^ We plot the redshift- 
dependence of the signal-to-noise ratios of the three methods in 
Fig. [5] We also plot the ratio of the signal-to-noise of the optimal 
polarisation method with the one of the standard method (solid) 



2 This will provide accurate results for multipoles 1 » 1, however, is a 
coarse approximation in the regime I ~ 1. 



3 The ratio of this neglected coupling to the template strength gets large for 
small 2 mlx . Our estimates are therefore less accurate in this regime. 
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Figure 5. Comparison of the signal-to-noise ratios versus the maximal red- 
shift 2 max of the galaxy survey. Top panel: Average signal-to-noise ra- 
tio of the optimal polarisation method (S/N)^^ (solid), of the optimal 
temperature-only method (S/N)f v (dashed) and signal-to-noise ratio of 
the standard method (S/7V) s t (dotted). Bottom panel: Ratio of the signal- 
to-noise of the optimal polarisation method with the one of the standard 
method (solid) and with the one of the optimal temperature-only method 
(dashed). We see that with polarisation data included, the signal-to-noise is 
significantly enhanced even for low redshifts. 



all secondary effects. However, in this work we have used a Gaus- 
sian approximation for the uncertainty in the signal template, which 
may not be optimal for effects on smaller scales, such as the RS ef- 
fect, the kinetic SZ effect or lensing. We leave the extension of our 
method to non-Gaussian noise models for future work. 

We calculated the achievable reduction in primordial noise for 
perfect (noiseless) data using the example of the ISW effect, and 
obtained a signal-to-noise ratio of up to 8.5. This corresponds to 
an enhancement of the signal-to-noise ratio by 16 per cent as com- 
pared to our optimal temperature-only method, independent of the 
depth of the LSS survey. In comparison to the standard method, 
the signal-to-noise ratio is enhanced by 23 per cent for a full-sky 
galaxy survey which goes out to a redshift of at least two. When 
using the SDSS main galaxy sample, which has a maximal redshift 
of about Zmax ~ 0.3, our signal-to-noise ratio is still enhanced by 
about 16 per cent as compared to the standard method. 

The variance reduction achieved with this method will sig- 
nificantly improve the detection of all kinds of secondary effects 
on the CMB, where a spatial template constructed from non-CMB 
data can be created. This stresses the importance of accurate mea- 
surements of primordial polarisation fluctuations even for non- 
primordial signal detection and analysis. The upcoming Planck 
Surveyor Mission, as well as more future experiments like Polar- 
BeaPJj or CMBPoQ will allow us to benefit from polarisation for 
the detection of secondary CMB signals in the way presented in 
this work. 



and with the one of the optimal temperature-only method (dashed). 
Note that the enhancement of the signal-to-noise ratio w.r.t. the op- 
timal temperature-only method is almost constant in redshift. This 
is quite clear from the fact that we have reduced the primordial 
noise with the polarisation data, and neither the primordial noise 
nor the reduction of the latter depend on redshift. Therefore, the 
reduction of the noise from including polarisation data is always 
the same, independent of how deep in redshift our survey goes, 
and the signal-to-noise ratio is already significantly enhanced for 
currently available surveys. For example, for a maximal redshift 
of Zmax ~ 0.3, which is the maximal redshift for the SDSS main 
galaxy sample, we have a better signal-to-noise by about 16 per 
cent as compared to the standard method. The additional enhance- 
ment for higher redshifts of our signal-to-noise ratio w.r.t. the stan- 
dard method comes from worki ng conditional on the g alaxy data, 
as we have described in detail in lFrommert et alj J2008l) . 



5 CONCLUSIONS 

The detection of secondary effects on the CMB remains a chal- 
lenge, because the amplitudes of these effects are much smaller 
than those of primordial CMB fluctuations. The techniques for de- 
tecting such secondary signals are all based on the existing cross- 
correlation between the LSS and the signal in question. However, in 
all of these studies, chance correlations of primordial CMB fluctua- 
tions with the LSS are the dominant source of noise in the analysis. 

We have presented a way of reducing the noise coming from 
primordial temperature fluctuations by simply subtracting the part 
of the temperature map which is known from the polarisation data. 
Effectively, only the unknown part of the temperature fluctuations 
then contributes to the variance of the signal estimate. 

As presented here, our method can be generically applied to 
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+ 



cr 



--( 



cr 



+ 1 



cr 



(A3) 



where we have completed the square in the second last step and 
used a generalised definition of the reduced temperature map and 
power spectrum, which we had introduced in eq. (122b . in the last 
step. Similarly, we can decompose the determinant of C(l): 



C(l) 



c AT ohs( jAE a 



cr 



Inserting eqs fA3l and <A4j into Q(af m — aj m , C(l)) allows us to 
write 



G(af m -a[ m ,C(l)) = G{af™ d - 
x G(af^ 



a lm:^l ) 
E T r AE 
a lm ' I 



(A4) 



In the case of the polarisation template E T being zero, this expres- 
sion reduces to the one in eq. l !23t . 



APPENDIX A: PROOF OF THE FACTORIZATION OF 
THE LIKELIHOOD 

We now explicitely prove the factorization of the likelihood in eq. 
d21t into a reduced temperature part and a polarisation part, as 
given in eq. {23}. We will do this for the more general case that 
we not only have a signal template T T for the temperature part, but 
also a non-zero template E T for the polarisation part. In this case, 
the covariance matrix is 

C(l) = c fi hMhs ° ( c AE obs J , (Al) 

instead of the simplified one given in eq. &201 . Here, A_E ob s is de- 
fined as Ai? ob s = £"obs — E T . The inverse of the covariance matrix 
is given by 

cor 1 = 1 



C fT obsC AE obs _ ^AT oba ,Afi obs j 



I 

AE obs _ r AT obs ,AE a 



_pAT ob „AE obB C ATo] 



(A2) 



We first rewrite the exponent of G(af m — a,] m , C(l)) in eq. J21b by 
inserting the inverse of C(l): 

,T T „B ob3 _ ^E T \ ^ A -l (T obs _ n T T n E obB _ n E T 



- 2 ( C r° bMb °/Cr^) Re 
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